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ABSTRACT 

The birth and death of the first generation of stars have important implications for the thermal 
state and chemical properties of the intergalactic medium (IGM) in the early universe. Sometime after 
recombination, the neutral, chemically pristine gas was reionized by ultraviolet photons emitted from 
the first stars, but also enriched with heavy elements when these stars ended their lives as energetic 
supernovae. Using the results from previous high-resolution cosmological simulations of early structure 
formation that include radiative transfer, we show that a significant volume fraction of the IGM can be 
metal-polluted, as well as ionized, by massive Population III stars formed in small-mass (~ 10 6 — 10 7 M Q ) 
halos early on. If most of the early generation stars die as pair-instability supernovae with energies up 
to ~ 10 53 ergs, the volume-averaged mean metallicity will quickly reach Z ~ lCU 4 /?© by a redshift 
of ~ 15 — 20, possibly causing a prompt transition to the formation of a stellar population that is 
dominated by low-mass stars. In this scenario, the early chemical enrichment history should closely trace 
the reionization history of the IGM, and the end of the Population III era is marked by the completion 
of reionization and pre-enrichment by z ~ 15. We conclude that, while the pre-enrichment may partially 
account for the "metallicity-floor" in high-redshift Lyman-a clouds, it does not significantly affect the 
elemental abundance in the intracluster medium. 

Subject headings: cosmology: theory — galaxies: formation — intergalactic medium — stars: formation 



1. INTRODUCTION 

Chemical elements heavier than lithium are thought to 
be produced exclusively through stellar nucleosynthesis. 
The primordial cosmic gas remains chemically pristine un- 
til the first supernova (SN) explosions expel metals that 
are produced in the precursor stars. High-redshift obser- 
vations of the Lyman-a forest (e.g., Songaila 2001; Pettini 
et al. 2003), damped Lyman-a systems (e.g., Prochaska 

2002) , and quasars (e.g., Freudling, Corbin, & Korista 

2003) all indicate that heavy elements are distributed in 
various cosmic environments by z ~ 5. While the ac- 
tual transport mechanism and the nature of the primary 
sources have not been determined, it is often suggested 
that the origin of these heavy elements may be attributed 
to the first generation of stars, the so-called Population III 
stars (e.g., Ostriker & Gnedin 1996; Madau, Ferrara, & 
Rees 2001; Thacker, Scannapieco, & Davis 2002; Qian, 
Sargent, & Wasserburg 2002). 

In the standard cosmological models based on cold dark 
matter (CDM), the first cosmological objects are predicted 
to form at redshifts greater than 20 (e.g., Couchman & 
Rees 1986; Tegmark et al. 1997; Abel et al. 1998; Yoshida 
et al. 2003a). Thus, the first heavy elements are likely 
to have been processed at such early epochs. Intriguingly, 
the recent measurement of a large Thomson optical depth 
by the Wilkinson Microwave Anisotropy Probe (WMAP) 
satellite (Kogut et al. 2003; Spergel et al. 2003) pro- 
vides evidence that the universe was reionized very early 
on, supporting the above notion that the first stars could 
form at z > 20 (e.g., Cen 2003; Haiman & Holder 2003; 
Wyithe & Loeb 2003a,b; Sokasian et al. 2003a). 



The relative abundances of heavy elements produced in 
the first stars are of great importance because observa- 
tions of the elemental abundance pattern of ultra metal- 
poor stars (e.g., Christlieb et al. 2002) can place strong 
constraints on the environments in which these stars were 
formed and possibly on the progenitor mass (Umeda & 
Nomoto 2003; Schneider et al. 2003a). Theoretical model- 
ing of the formation of the first stars (Abel, Bryan, & Nor- 
man 2002; Bromm, Coppi, & Larson 2002; Omukai & Palla 
2003) consistently indicate that they were rather massive, 
with characteristic mass being possibly several hundred 
solar masses. If the first stars are indeed as massive as 
~ 200M©, they end their lives as energetic SNe via the 
pair-instability mechanism (e.g., Barkat, Rakavy, & Sack 
1967; Bond, Arnett, & Carr 1984; Fryer, Woosley, & Heger 
2001; Heger & Woosley 2002; Woosley, Heger, & Weaver 

2002) , releasing a total energy of up to ~ 10 53 ergs. Such 
energetic explosions in the early universe are violently de- 
structive: they expel the ambient gas out of the gravi- 
tational potential well of small-mass dark matter halos, 
causing an almost complete evacuation (Bromm, Yoshida, 
& Hernquist 2003, hereafter Paper I; Wada & Venkatesan 

2003) . 

Since the massive stars process a substantial fraction of 
their mass into heavy elements, early SN explosions may 
provide an efficient mechanism to pollute the surround- 
ing intergalactic medium (IGM). Models of stellar evolu- 
tion (e.g., Heger & Woosley 2002) predict that, in massive 
zero-metallicity stars, up to ~ 40 — 50% of the total stel- 
lar mass is processed into heavy elements in the core and, 
in the pair-instability case, all of it is finally ejected dur- 
ing the explosion without leaving any compact remnant 
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behind. Blast waves triggered by the SN explosions will 
then drive supersonic flows of gas that are enriched with 
the processed heavy elements. Consequently, because of 
the short lifetimes of such massive stars, the emergence of 
the early generation stars should be almost immediately 
followed by metal enrichment of the IGM; i.e., prompt 
metal enrichment could be achieved efficiently by Popula- 
tion III SNe in the early universe. 

The transport of metals by SN-driven winds has been 
studied extensively in the context of conventional galaxy 
formation (e.g., Larson 1974; Dekel & Silk 1986; Vader 
1986; Mac Low & Ferrara 1999; Aguirre et al. 2001a,b; 
Madau ct al. 2001; Scannapicco, Thacker, & Broadhurst 
2001; Mori, Ferrara, & Madau 2002; Furlanetto & Loeb 
2003; Springel & Hcrnquist 2003a, b). In fact, gas flows 
induced by SN explosions appear to be the most plausible 
way of distributing metals, which are observed to exist un- 
der a broad range of environments. Interestingly, an array 
of observations indicate that the metals were produced and 
dispersed at early stages of structure formation. Spectra 
of distant quasars consistently show that there are at least 
some metals in the diffuse IGM at z ~ 2 — 3 (e.g., Schaye 
et al. 2000, 2003). Qian & Wasserburg (2001, 2002) ar- 
gue that the observed abundance patterns of metal-poor 
stars can be explained by the "prompt-inventory" of heavy 
elements from very massive stars. Loewenstein (2001) 
points out that a contribution from very massive metal- 
free stars may account for the observed abundance anoma- 
lies of oxygen, silicon, and iron in the intracluster medium 
(ICM). More recently, Frcudling et al. (2003) discovered 
a strong Fe II emission feature from a z ~ 6 quasar, 
providing strong evidence that iron is produced in some 
stars/galaxies already at z J> 10, which is expected if black 
hole growth is regulated by star formation (e.g., Di Matteo 
et al. 2003a, b). Together with the implication of an early 
reionization derived from the WMAP results, these obser- 
vations indicate that Population III stars may have been 
an important contributor of heavy elements in the diffuse 
IGM/ICM as well as the source of radiation responsible 
for the early reionization of the universe. The association 
of the IGM metal enrichment with the UV-photon produc- 
tion from the first stars was previously discussed by Oh ct 
al. (2001), and more recently by Stiavelli, Fall, & Pana- 
gia (2003). In these studies, however, the analysis was 
not based on a detailed modeling of the high-redshift star 
formation rate, and the precise epochs of Population III 
pre-enrichment and reionization therefore remained some- 
what uncertain. 

In the present paper, we revisit the cosmic metal pro- 
duction in the context of early structure formation. We ex- 
plore the possibility that the first stars are responsible for 
reionization of the universe at z ^ 15, while also causing a 
significant, prompt chemical enrichment of the IGM. We 
argue that, if the early reionization was the result of very 
massive Population III stars, the mean metallicity should 
reach a critical level at which the efficiency of gas cooling 
is greatly enhanced, thus changing the overall mode of star 
formation (e.g., Mackey, Bromm, & Hcrnquist 2003). 

The hydrodynamic transport of heavy elements induced 
by SN explosions is studied numerically in Paper I, us- 
ing high-resolution, three-dimensional cosmological simu- 
lations. We show that a large amount of the gas is ex- 



pelled out of the shallow gravitational potential of the host 
"minihalos" . The simulations in Paper I further demon- 
strate that, if the explosion kinetic energy is as high as 
~ 10 53 ergs, nearly 90% of the pre-enriched gas in the 
immediate vicinity of the explosion site eventually escapes 
out of the host halo, reaching up to a distance of ~ 0.5 kpc 
in a few million years. Since multiple SNe are expected 
to be occurring at z ^ 15 in distant star- forming regions 
within a cosmological volume, the total amount of metals 
produced and distributed into the IGM may well exceed 
the critical level required to affect the thermal and chem- 
ical properties of the IGM. The transition in stellar pop- 
ulations, from Population III to Population II, may then 
be naturally caused early on. 

Throughout the present paper, except in section ?? 
where we consider a variant CDM model, we work with 
a flat A-dominated Cold Dark Matter universe with mat- 
ter density £l m = 0.3, cosmological constant Ha = 0.7 
and the Hubble constant at the present time h = 0.7 in 
units of 100km s^Mpc -1 . We set the baryon density to 
£lb = 0.04 and the initial density fluctuations are normal- 
ized to og = 0.9. These parameters are consistent with 
those employed in the numerical simulations we refer to in 
the following sections. 

2. FEEDBACK FROM EARLY GENERATION STARS 

In this section we summarize various feedback effects 
from the first stars, and derive useful order-of-magnitude 
estimates for quantities needed to study these effects. The 
ultimate fate of a metal-free star depends critically on its 
mass. Theoretical stellar models predict the following clas- 
sification (Heger & Woosley 2002; Heger et al. 2003): (1) 
8 < M» < 25M Q : these stars explode as core-collapse SNe 
and leave a neutron star behind, (2) 25 < M* < 40M Q : 
these explode as faint Type II SNe and leave black holes 
behind after fall-back of most of the envelope (see Umeda 
& Nomoto 2002, 2003), (3) 40 < M* < 140M Q : these do 
not explode as SNe and directly collapse into black holes, 
(4) 100 < < 140M Q : these experience a pulsational 
instability and eject their outer envelope, again leaving a 
black hole behind, (5) 140 < M* < 26OM : these explode 
as pair-instability SNe, causing complete disruption, and 
(6) M* > 26OM : these collapse, in the absense of rota- 
tion, directly into massive black holes. 

In summary, only the pair-instability case can substan- 
tially contribute to the metal enrichment of the IGM, 
whereas remnant black holes (or neutron stars) are left 
in the other cases, which we shall discuss in § 2.5. While 
early generation stars may be born with a broad range of 
masses (e.g., Bromm, Coppi, & Larson 1999; Nakamura 
& Umemura 2001, 2002; Omukai & Yoshii 2003), very 
massive stars are dominant in causing various feedback 
effects to the IGM in the early universe such as ioniza- 
tion, heating, and chemical enrichment. Thus, through- 
out the present paper, we primarily consider very massive 
(~ 100 — 3OOM ) metal-free stars as Population III stars. 

2.1. Population III Stars as Reionization Sources 

We first discuss radiative feedback. The efficiency of 
ionizing photon production from massive Population III 
stars is obtained from detailed calculations (e.g., Bromm, 
Kudritzki, & Loeb 2001b; Schaercr 2002, 2003). As our 
fiducial case, we consider a metal-free star with mass 
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M* = 200M Q . Such a star emits about 3 x 10 50 ioniz- 
ing photons per second, and has a main-sequence lifetime 
of ~ 2 million years. 

Consider a population of such massive stars formed in a 
cubic comoving-Mpc volume by z. We introduce a quan- 
tity 

•M = (i) 

where iV p h(z) is the cumulative number of ionizing photons 
emitted from the stars, and Nn = 4 x 10 66 is the number 
of hydrogen atoms in the volume. This 'number of pho- 
tons per atom' q(z) serves as an indicator of the global 
ionization fraction as a function of z. Clearly, in order to 
completely ionize the IGM, at least the same number of 
UV photons as the number of hydrogen atoms must be 
emitted within the volume. It follows from these numbers 
that more than a few hundred very massive stars must 
be formed in a comoving-Mpc volume to cause complete 
reionization. 

2.2. Population III Stars as Sources of Heavy Elements 

Massive metal-free stars are not only luminous ionizing 
sources, but also efficient factories of heavy elements. For 
stars with masses ~ 100 — 3OOAf , typically up to ~ 50% 
of the mass is processed into heavy elements in the core, 
consisting mostly of carbon, oxygen, silicon and iron. We 
denote this mass fraction with fz, and set it to be 0.5 as 
our fiducial value for massive stars in the remainder of the 
paper. We further introduce a parameter 



s N~, 
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where N n is the total number of nuclei in the heavy ele- 
ments processed in the core, and 7V 7 is the total number 
of ionizing photons emitted by the star during its lifetime. 
Namely, with m p denoting the proton mass, m p x N n is the 
total mass of metals produced in the star. This quantity 
gives the ratio of the stellar metal to photon production 
which is, in principle, uniquely determined by the struc- 
ture and evolution of the respective star. For a metal-free 
star with mass 200Af Q that is destined to die in a pair- 
instability supernova, the actual value is approximately 
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fz x 2OOM /m p 
2 x 10 64 



10 



r,9 



2 x 10 64 



5 x 10~ 6 , (3) 



f° r fz — 0.5. It is straightforward to compute the value for 
£ for less massive stars or for stars with finite metallicity 
from stellar evolution models. 

Intriguingly, this parameter Cpisn for the very massive 
stars that are destined to explode as pair-instability SNe 
can be expressed in terms of fundamental physical con- 
stants as follows. We estimate the total number of ion- 
izing photons produced during the lifetime of a Popula- 
tion III star with mass M* as iV 7 ~ 0.007M*c 2 /ei on - Here, 
ei on ~ e 2 /ao — a 2 m e c 2 is the binding energy of an electron 
in an H atom, with ao being the Bohr radius, a ~ 1/137 
the fine-structure constant, and m e the electron rest mass. 
Assuming the limiting case of a 100% nuclear burning effi- 
ciency, we can write for the total number of nucleons pro- 
cessed into heavy elements: N n ~ M*/m p . We now have, 



for very massive stars that die as pair-instability SNe, 

CplSN = 4^L~4xnr 6 , (4) 



where a 2 n ~ 0.007 describes the strength of the pion- 
nucleon interaction. The production of ionizing photons 
and metals in a massive star is thus determined by the 
ratio of the atomic and nuclear binding energies. 

2.3. Population III Supernova Explosions 

The processed metals can be distributed into the IGM 
only if some transport mechanism operates, otherwise the 
metals would be confined within a small region and the 
diffuse IGM would be largely unaffected. In this section, 
we first describe the hydrodynamics of a single Population 
III SN explosion in the early universe and then show that 
the gas and the heavy elements can be transported over a 
large distance. 

Assuming spherical symmetry of the system, we apply 
the usual Sedov- Taylor self-similar analysis for a given ex- 
plosion energy and density of the ambient medium. We 
first need to specify the density profile of the surround- 
ing IGM in order to fully describe the evolution of the SN 
blast wave. For a density profile approximated by a sim- 
ple power-law, p ex r~ w , the expansion of a blast wave is 
described by the self-similar solution: 



+2\ 1/(5-™) 



(5) 



where £ is a constant determined by the energy integral, 
E is the explosion energy, n the particle number density of 
the surrounding medium, and t the time elapsed since the 
explosion (Sedov 1959; see also Ostriker & McKee 1988). 
In reality, the gas density profile can be more complex be- 
cause photo-evaporation due to radiation from the central 
star is very effective in small halos (T v ; r <^ 10 4 K). Large 
pressure gradients in the ionized region drive an outward 
gas flow before the explosion is triggered (see Paper I). 
Hence, although initially the gas density is high around 
the star, it is expected to be substantially reduced by the 
time when the star explodes, and the density profile will 
be much shallower except in the outer edge where a su- 
personic Shockwave is propagating. Taking the effective 
photo-evaporation into account, we assume the density of 
the ambient medium around the explosion site is approxi- 
mately constant at the mean IGM density at z, 



hh.igm ^ 0.0014 
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20 



(6) 



where fib is the baryon density and h is the Hubble con- 
stant in units of 100 km s _1 Mpc^ 1 . Note that the shock 
radius in equation (J5J scales only weakly with the density 
and the explosion energy unless the density profile is very 
steep. The three-dimensional simulations in Paper I show 
that the evolution of the shock radius is well described by 
equation 10 in a piecewise manner for regions having den- 
sity profiles with different slopes. For a constant density 
medium, the Sedov- Taylor self-similar solution reduces to 



R s ~ 200 
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pc, 
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where R s is the shock radius, -E51 the total SN energy in 
units of 10 51 ergs, n the particle number density of the 
ambient medium in units of 1 cm~ 3 , and t? the elapsed 
time since the explosion in units of 10 7 years. For the 
explosion energy, Heger & Woosley (2002) calculate the 
energy released by a Population III pair-instability SN to 
be in the range ~ 10 51 — 10 53 ergs. 

The above solution is valid if the gas remains adiabatic, 
but breaks down when the shock becomes radiative. The 
hot gas in the interior of high redshift (z > 10) SN rem- 
nants (SNRs) cools partly via inverse Compton cooling. 
The characteristic cooling time for this process is given by 

*Comp =s 10 7 (^p) y- ( 8 ) 

Note the strong dependence on redshift. At z = 20, the 
SNR will remain approximately adiabatic (hence equation 
jZj is valid) until ~ 10 7 years after the explosion. How- 
ever, this provides only a maximum timescale. In the post- 
shock region, the gas density and temperature are so high 
that hydrogen and helium atomic line cooling is the dom- 
inant cooling process. Thus, a thin dense shell could be 
formed faster than the above timescale. The exact time 
when the post-shock region becomes radiative depends on 
the total explosion energy and the density of the ambient 
medium in a complicated manner. The shell forms when 
the local cooling time i C ooi(?~) at radius r becomes shorter 
than the expansion timescale i? s /i? s . We use the ana- 
lytic Sedov- Taylor solution for the temperature T(r) and 
density p(r) of the gas to compute the local gas cooling 
rate A(p(r), T(r)), assuming an explosion energy E — 10 53 
ergs and the gas density n = 0.0014 cm~ 3 (~ the mean 
IGM density at z ~ 20). Recent numerical simulations by 
Whalen, Abel, & Norman (2003) indeed show that photo- 
evaporation by the central star reduces the final gas den- 
sity within the host halo to very small values, comparable 
to the mean IGM density. In this case, the shell forms at 
a time ~ 3 x 10 6 yr. Then the shell radius (the "stalling- 
radius") is ~ 1.2 kpc, and it expands only slowly during 
the later snow-plow stage as R s oc t 1 ' 4 (see Figure 1 in 
Paper I). This is perhaps an optimistic estimate because 
the gas evacuation efficiency by radiation from the stars 
depends on the host halo mass and the gas density pro- 
file. For a ten times higher density, n = 0.014 cm -3 , the 
stalling-radius is calculated to be only 40% smaller, ~ 0.7 
kpc, which is more consistent with the result of our simu- 
lation presented in Paper I. Overall, because of the weak 
dependence of the self-similar solution (equation [7]) on 
the density, the above estimates will not be significantly 
affected by (possible) variations in the gas density around 
the SN sites. 

The momentum- driven slow expansion is still impor- 
tant over a cosmological timescale (Bertschinger 1983), 
with the radius doubling every twenty folds in time. 
Thus, if mixing of the original ejecta with the ambient 
medium is efficient, the expelled metals can travel over 
a large distance (~ a few kpc) within a Hubble time. 
Three-dimensional numerical simulations of Nakasato & 
Shigeyama (2000) show that the SN ejecta are fairly well 
mixed within the remnant, with the enriched volume frac- 
tion being close to unity. Wada & Venkatesan (2003) show 
that metal-dispersal is very efficient even in dense, disk- 



shaped gas clouds. Therefore, the processed metals are 
likely to be transported as far as the shell radius. 

2.4. Global Quantities 

Let us calculate the evolution of the volume-averaged 
mean metallicity of the IGM. We allow the uncertainty 
associated with the characteristic mass, or initial mass 
function (IMF), of Population III stars by assuming that 
a fraction, /pisn, of the Population III stars have masses 
in the pair-instability range. Suppose that n*(z) stars are 
formed by z in a unit comoving volume. The total mass 
of metals produced within the volume is then 

M z =/pisNXfi,(i)/ z M«. (9) 

We further normalize this quantity to the mass of the IGM 
within the volume, to obtain the mean metallicity: 

2 _ M z _ njjz) x /pi SN x fzM* 
Migm m p x iV H 

~ /pisn x q(z) x Cpisn- (10) 

Recall that Cpisn is approximately a constant for mas- 
sive metal-free stars that explode as pair-instability SNe 
(equation [3]), and q(z) traces the IGM ionization fraction. 
The above equation simply shows that the total quantity 
of heavy elements ejected into the IGM is proportional 
(weighted by the fraction /pisn) to the number of UV 
photons emitted by z. 

We now consider the simplest case where all of the Pop- 
ulation III stars have a mass of 200Mq and explode as 
pair-instability SN, i.e., /pisn=T- We further assume that 
9reion ~ 2 marks the epoch of reionization (see, e.g., fig. 8 
of Sokasian et al. 2003b). Setting Cpisn = 5x 10 -6 , we 
obtain 

%,„«,■> = 10_5 ~ 5 x !0~ 4 Z©, (11) 

which coincides with the critical metallicity necessary for 
the formation of low-mass gas clumps (Bromm et al. 
2001a; see also Mackey et al. 2003). Note also that very 
massive stars become progressively unstable for increas- 
ing metallicity (Baraffe, Heger, & Woosley 2001). There- 
fore, for /pisn = 1) the epoch when the transition in the 
mode of star formation occurs should closely match the 
first reionization epoch. Note, however, that we have ne- 
glected recombination in the above argument. If Popula- 
tion III star formation occurs rather gradually, taking a 
considerably longer time than a recombination time, most 
of the UV photons are simply consumed by recombina- 
tion without increasing the global ionization fraction of 
the IGM. Approximating the IGM density with the cosmic 
mean density, we estimate the recombination timescale for 
hydrogen: 

* rcc «i0 8 yr(i±£) 3 - (12) 

If the increase in the number of UV photons is sufficiently 
rapid, compared to t rcc , then the ionization fraction of the 
IGM should monotonically increase. This is indeed the 
case for the model we describe below. 

An important question is when and how rapidly did 
the early generation stars emerge? We address this is- 
sue using the formation rate of Population III stars de- 
rived by Yoshida et al. (2003a). Their result is based 
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on a semi-analytic model coupled with the outputs of a 
large cosmological A-body simulation. In their model, 
star-forming regions are identified as halos in which the 
gas has been able to cool via molecular hydrogen cool- 
ing, and the usual 'one- star- per- halo' assumption is made 
such that only a single star is formed per star-forming 
gas cloud. The model also includes the negative feedback 
owing to photo-dissociation of molecular hydrogen by pho- 
tons in the Lyman- Werner band. Hence, the model pro- 
vides a conservative estimate for the early star-formation 
rate (SFR). It is worth noting that the resulting SFR was 
found to be in good agreement with that obtained from 
the hydrodynamic simulations of Sokasian et al. (2003a) 
which include radiative transfer as well as various feed- 
back processes. We revised the empirical fit to the SFR in 
Yoshida et al. (2003a) and found that it is well-described 
by a particular form proposed by Hcrnquist & Springel 
(2003): 



P*{z) = P 



X" 



*.o 



l + 5(x-l)9exp(/? X r ) : 



where 



X(z) 



H(z] 



2/3 



(13) 



(14) 



and, in this case, p = 5, q = 1, r = 1.85, a = 0.01, 
/? = 0.04, p\ >0 = 3.0 x 1O- 8 M yr- 1 Mpc" 3 . 

Figure 1 shows the Population III SFR. The flattening 
at z < 25 reflects the suppression of the gas cooling due 
to global radiative feedback, showing the self-regulating 
nature of early star formation (e.g., Ricotti, Gncdin, & 
Shull 2002; Wise & Abel 2003; Yoshida et al. 2003a). We 
compare the derived SFR with the result of the hydrody- 
namic simulations of Sokasian et al. (2003a) . To this end, 
we have averaged over all the nine models considered in 
Sokasian et al. (2003a) which differ in the star-formation 
efficiency as well as the strengths of the feedback effects. 
Overall, the agreement between the semi-analytic predic- 
tion and the simulation results is quite good. As Sokasian 
et al. (2003a) note, the SFR derived from their simula- 
tions may suffer from finite boxsize effects at high red- 
shifts (z > 25), which accounts for the small discrepancy 
between the analytic model and the simulation result seen 
in Figure 1 (see also Yoshida et al. 2003b; Barkana & Loeb 
2003) 

In the middle panel of Figure 1 , we plot the cumulative 
number of ionizing photons per baryon q(z). We compute 
the number of ionizing photons assuming a single massive 
(M* = 200M Q ) star emits ionizing photons at a rate of 
3 x 10 50 s _1 during its main-sequence lifetime (Bromm et 
al. 2001b). The dashed line in the middle panel indicates 
q = 1, which is achieved by z — 17. We note that the 
total number of UV photons increases by more than an 
order of magnitude between z — 22 and z = 17. The 
corresponding time interval is ~ 50 Myr, whereas the re- 
combination timescale at z ~ 20 is ~ 100 Myr. Therefore, 
to lowest order, the photon consumption due to hydrogen 
recombination can be neglected. In addition, since the 
photon escape fraction from minihalos is expected to be 
of order unity (Whalen et al. 2003), and the clumping of 
the IGM at these early epochs is small (Haiman, Abel & 
Madau 1999), q ~ a few should suffice to reionize a sub- 
stantial fraction of the IGM (e.g., Sokasian et al. 2003b). 



It is worth mentioning that, in our model, the time when 
q = 1 coincides with the epoch of rcionization suggested by 
the WMAP data, z = 17 ± 4 (assuming an instantaneous 
reionization, sec Kogut et al. 2003). At the same time, as 
shown in the bottom panel in Figure 1, the global mean 
metallicity reaches the metallicity floor of a few x \Q~ a Zq. 
We emphasize that this global mean is merely an approxi- 
mate estimate because, locally around the SN sites where 
the formation of second generation stars is likely to take 
place, the gas metallicity could be much higher than this 
value. 
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Fig. 1. — The global star formation rate density (top), the cu- 
mulative number of photons per hydrogen atom (middle), and the 
mean IGM metallicity Z (bottom), as functions of redshift. In the 
top panel we compare our semi-analytic prediction (solid line) with 
the result of numerical simulations of Sokasian et al. (2003a) (filled 
circles). The dashed line in the middle panel indicates q = 1. The 
critical metallicity 10~ 3,5 Zq is also indicated by the dashed line in 
the bottom panel. We plot three cases with mixing efficiency of the 
heavy elements with the diffuse IGM, 100% (solid line), 50% (dotted 
line), and 5% (dot-dashed line). 



2.5. Remnant Black Holes 

In the previous sections, we have exclusively considered 
very massive stars that explode as pair- instability SNc, i.e., 
those with masses 140 — 26OM . Stars with masses outside 
this range are likely to leave remnant black holes behind 
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(Fryer et al. 2001). In particular, stars even more mas- 
sive than 260Af Q completely collapse to black holes, which 
could be the seeds of supermassive black holes (Madau & 
Rees 2001; Islam, Taylor, & Silk 2003; Volonteri, Haardt, 
& Madau 2003). 

It is straightforward to compute the evolution of the 
mass density locked up in black holes using our model Pop- 
ulation III SFR. Let us adopt the simplest model where 
all the Population III stars are extremely massive and be- 
come black holes. Although this is an implausible situation 
because then the IGM would remain chemically pristine 
forever and the "star-formation conundrum" would arise 
(see Schneider et al. 2002), it is illustrative to study the 
evolution of the black hole mass density in this scenario. 

We consider three cases with M« = 300, 600, IOOOMq, 
and the resulting mass density is plotted in Figure^ Note 
that our model for the Population III star formation is 
based on the 'one-star-per-halo' assumption and thus the 
star formation rate density scales approximately with M * . 
Figure|21shows that the mass density at z ~ 15 is compara- 
ble to the present-day mass density in black holes in the lo- 
cal universe (e.g., Merritt & Ferrarese 2001) for high mass 
values. Note that these Population III blackholes are likely 
to have sunk to the galactic centers by the present epoch 
due to dynamical friction (e.g., Madau & Rees 2001). 

Obviously, this extreme top-heavy mode of Population 
III star formation should be stopped at some point. Oth- 
erwise too much baryonic mass would be converted into 
black holes locally by z = 0. 




redshift 

Fig. 2. — The total mass of baryons converted to stars per unit 
comoving volume, for M* = 300 (solid line), 600 (dashed line), and 
1000 Mq (dot-dashed line). The horizontal dotted line indicates 
the present-day mass density in black holes in the local universe 
estimated by Merritt & Ferrarese (2001). 

3. IMPLICATIONS 
3.1. The End of the Population III Era 

One may naively expect that Population III star forma- 
tion could continue down to moderate redshifts, z ~ 5 — 10. 
We argue, based on the results of our model, that either 
the completion of reionization or prompt enrichment to a 
level of Z ~ 10~ 4 Z Q by z ~ 15 should mark the end of 
the Population III era. 

When a large fraction of the IGM is reionized, effec- 



tive photo-evaporation of minihalos by internal and ex- 
ternal sources occurs (e.g., Shapiro, Iliev, & Raga 2003), 
and virtually no gas can be trapped within low mass ha- 
los. Furthermore, SN explosions cause nearly complete 
evacuation of the host halos (see Paper I). Therefore, re- 
gardless of the efficiency of the proposed positive feedback 
on molecular hydrogen formation (e.g., Ferrara 1998; Oh 
2001), it should take some time for the heated and ex- 
pelled gas to cool and condense. Subsequent star forma- 
tion can take place only after the gas falls back into the 
halos again. Accordingly, there should be an appreciable 
time gap between the first star formation and the onset of 
the second-generation star formation. We conservatively 
estimate this 'recuperation time' (Sokasian et al. 2003a) 
to be of the order of the dynamical time scale for the halos. 
The dynamical time scale is given by 

* dy * = \^fp ' (15) 

where G is the gravitational constant and p the total mass 
density of a halo. Assuming that the mass density of a halo 
(including dark matter) is 200 times the cosmic mean den- 
sity, we compute the dynamical time for halos collapsing 
at z = 15 to be ~ 50 Myr. Since the cosmic time differ- 
ence between z — 15 — 18 (or z — 13 — 15) is ~ 60 Myr 
for the standard ACDM model, star formation is likely 
to be quenched for at least a duration of Az ~ a few in 
redshift. It might seem that halos with virial tempera- 
ture greater than 10 4 K are unaffected by photo-heating 
due to reionization (e.g., Oh & Haiman 2002). We point 
out that, in the hierarchical universe, a^the massive halos 
used to be minihalos at earlier epochs. Hence, according to 
our model of Population III star-formation, the gas within 
them is likely to have been already expelled either by ra- 
diation from the Population III stars or by SN explosions. 
The situation would be slightly different in halos that are 
rapidly growing in mass because such halos may be able 
to trap the ionized hot gas. However, SN explosions with 
energies ~ 10 53 ergs are still able to cause an outward 
gas flow even in the large halos. If the feedback from the 
first massive star on its vicinity occurs by z=15, as in our 
model, then the formation of the second-generation stars 
in the same region (i.e., protogalaxy) will begin at z < 13. 
Further, very massive stars are already pulsationally un- 
stable when the metallicity is ~ 10 _4 Zq, even if they are 
formed (Baraffe et al. 2001). The minimum metallic- 
ity floor is established in all the star-forming regions by 
z ~ 15 if the first stars explode as pair-instability SNe, 
thus causing a 'global' regulation of the star formation ac- 
tivity lasting until z ~ 13. Note that significant cosmic 
variance may be induced by large-scale (^> IMpc) density 
fluctuations. The cosmic scatter in the number density 
of the star-forming clouds should then cause a substantial 
spread in the epoch of self-termination (e.g. Scannapieco, 
Schneider & Ferrara 2003). 

An alternative scenario for second-generation star for- 
mation is self-propagating star formation in the thin, dense 
shell formed by early SN explosions (e.g., Ostriker & Cowie 
1981; Wandel 1985; McCray & Kafatos 1987). Such phe- 
nomena are indeed observed in the local universe (e.g., 
Dopita et al. 1985; Efremov & Elmegreen 1998) and may 
have taken place in the early universe as well (Mackey et 
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al. 2003; Salvaterra, Ferrara, & Schneider 2003). For the 
particular case we simulated in Paper I, we found that 
no gas clumps are formed in the dense shell which are 
Jeans-unstable. This may, however, be because of the still 
insufficient mass resolution of the simulation. In order to 
properly address the question of the shell fragmentation 
and the formation of proto-stellar clouds, simulations with 
an extremely high mass resolution are needed. Also, more 
complex configurations, such as collisions between neigh- 
boring blast waves originating in more massive systems, 
may result in the formation of gravitationally bound frag- 
ments. We pursue this possibility further in a forthcoming 
work (Bromm et al., in preparation). 

3.2. Thomson Optical Depth and Mean Metallicity 

In § 2, we have shown that the production of ionizing 
photons and heavy elements by massive Population III 
stars is described by fundamental constants in the case 
where these stars explode as pair-instability supernovae. 
Here, we derive an approximate but concise relationship 
between the evolution of the mean metallicity and the to- 
tal Thomson optical depth, r e , as follows. 

The differential optical depth to Thomson scattering 
corresponding to the redshift interval dz is 



dr e = (It n e (z) c— dz, 
dz 



(16) 



where or = 6.65 x 10~ 25 cm 2 is the Thomson scattering 
cross section, and n e (z) the mean electron number density 
at z. The mean electron number density is given by 



n e (z) = n e (0)(l + z) 3 x e (z), 



(17) 



where x e is the ionization fraction at z, and n e (0) the 
present-day number density of electrons. We obtain the 
contribution from an interval between z\ and z-i (z\ < Z2), 



dr e \ Zl ^ Z2 =ct t n e (0)x e (z) — 
(l + zf 



(18) 



(tt A + n m (l + z)3)l/2 



dz 



■ ct t n e (0) x e (z) 



c / 1 + z 



Ho V fl m 



1/2 



We have assumed a flat ACDM model with H(z) = 
Hq[Q\ + £l m (l + z) 3 ] 1 / 2 and the last approximation holds 
at high redshifts, certainly for z > 10. 

Neglecting recombination, we can approximate, to the 
lowest order, the ionization fraction as 

Xe(z) oc q(z). (19) 

From equation I jlOjl we obtain for pair-instability SNe, 

Z(z) 



x e {z) 



Cpisn/pisn 



(20) 



Absorbing all the constant factors into a single quantity 
T, including /pisn: we express the concise relation as 



dr„ 



T / Z(z)(l + z) l/2 dz. 



(21) 



Thus, to the lowest order, the evolution of the ionization 
fraction should be accompanied by the increasing metal- 
licity following equation 121|) . Integrating equation l|21|) to 
z = Zchom, when Z = 5 x 10 _4 Z Q , we obtain 



Ar c ,p piii ~ 0.06. 



(22) 



This is the contribution from early Population III stars 
expected in our model. 

3.3. Metallicity in the Lyman-a Forest at z ~ 3 — 5 

Metals produced by the early Population III stars should 
exist ubiquitously in various systems at low redshift. 
Miralda-Escude & Rees (1997) argue that the observed 
metallicity in the Lyman-a forest at z ~ 3 is consis- 
tent with a simple picture in which the metals originate 
from an early generation of stars. The average metallic- 
ity Z = lO -2 iv0 in the Lyman-a forest at z ~ 3 is still 
much higher than the mean metallicity achieved at the 
end of the Population III era, z ~ 15, whereas Songaila 
(2001) argues that a minimum universal metallicity of 
about 3.5 x V^T^Zq is found in z ~ 5 Lyman-a clouds. 
Remarkably, this is very close to Z ~ 3 x 10 _4 Z© (Fig- 
ure 1) achieved by the early massive Population III stars. 
Since the Population III star formation is self-regulating 
owing to the global radiative feedback and metal enrich- 
ment, this value is expected to be approximately the max- 
imum which early Population III stars alone can achieve. 
Songaila (2001) estimates the mass density of C IV be- 
tween z = 1.5 and z — 5 to be ^Civ.igm = (2.5 — 7) x 10~ 8 . 
Assuming C IV is the dominant ionization stage with ion- 
ization fraction of 0.5 (Songaila 1997), the mass density 
in carbon is ~ (5 — 14) x 10~ 8 . For a very massive zero- 
metalicity star, the carbon yield is ~ 4M© and nearly con- 
stant for the relevant core mass (Heger & Woosley 2002). 
From this we obtain the mass density in carbon (more 
precisely, 12 C) for our model with /pisn = 1, 



n 



Civ.PopIII 



10' 



(23) 



at z = 15. Although smaller than the value that Songaila 
(2001) obtains for the diffuse IGM at z ~ 5, this is not neg- 
ligible. Indeed, a fraction of the carbon, up to 20%, in the 
IGM may originate from the early generation Population 
III stars. 

At z ~ 3, the bulk of the metals found in the Lyman- 
a forest may also come from ordinary stellar populations 
formed in galaxies at lower redshift, z < 5 (e.g., Adel- 
berger et al. 2003). Interestingly, the metallicity evolution 
computed by Hernquist & Springel (2003) based on their 
analytic model of star-formation in galaxies predicts that 
the average metallicity reaches ~ 10~ 2 Zq at z ~ 3 — 5, in 
good agreement with the observed value in the Lyman-a 
forest. 

3.4. Enrichment of the Intracluster Medium 

X-ray observations of galaxy clusters reveal that the 
ICM is enriched to a significant fraction of the solar metal- 
licity (e.g., Mushotzky et al. 1996; White 2000). In 
fact, a large fraction (~ 75%) of a typical cluster's en- 
dowment with heavy elements exists in the ICM rather 
than in cluster galaxies (Renzini 1997). Clusters retain 
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all the metals produced in them, and the determination 
of the abundances for several elements in the ICM is 
relatively straightforward, without involving complicated 
modeling. Thus, the elemental abundances in the ICM 
offer a clean template against which any theoretical mod- 
els can be tested. Loewenstein (2001) suggests that pre- 
enrichment by a generation of massive Population III stars 
may account for the abundance anomalies of oxygen, sili- 
con and iron, observed in the ICM. The recent analysis of 
the relative overabundance of Si with respect to Fe in the 
ICM also indicates the possible contribution from massive 
Population III stars. Baumgartner et al. (2003) argue 
that the observed relative overabundance of Si can be ex- 
plained by Population III progenitors with a 70 M Q He 
core. While including the contribution from massive Pop- 
ulation III stars is attractive, the required abundance ap- 
pears very large. Baumgartner et al. (2003) estimate that 
the total mass in Population III stars should be at least 
^Popin ~ 10~ 5 , in units of the cosmic density parameter, 
whereas our model predicts only 

fipopin ~5x 10" 7 , (24) 

at z = 15. Baumgartner et al. (2003) argue that pri- 
mordial star formation is enhanced in overdense regions 
which later become clusters. However, while this might 
seem plausible, the strong self-regulation of early star for- 
mation primarily tends to shift star formation to earlier 
epochs in overdense regions, rather than locally enhancing 
the total amount of formed stars. We therefore conclude 
that metals originating from the earliest generation stars 
cannot be responsible for the observed abundance anoma- 
lies in the ICM. 

3.5. Dust Formation and H2 Production 

We have computed the mean metallicity of the IGM as- 
suming implicitly that all the metals produced and ejected 
by SN explosions will remain in the gas phase. This is not 
strictly true; a fraction of the metals will be in the form 
of dust and hence do not directly affect the gas chemi- 
cal properties (see however the effect of promoted molec- 
ular hydrogen formation discussed later in this section). 
Schneider, Ferrara, & Salvaterra (2003a) and Nozawa et 
al. (2003) calculate the dust depletion factor for pair- 
instability SNe to be up to 30-40%. Nozawa et al. (2003) 
also conclude that the dust formation rate and the relative 
abundances of the dust grains of various compounds de- 
pend sensitively on the mass of the progenitor and on the 
mixing efficiency in the ejecta. Another complexity is that 
the produced dust grains are destroyed by thermal sputter- 
ing with ions and electrons, and by grain-grain collisions 
in post-shock regions (Spitzer 1978). Most important, the 
reverse shock formed in the early stages of a SN explo- 
sion can destroy a significant fraction of the dust grains 
produced in the ejecta (McKee 1989). Unfortunately, the 
efficiency of dust destruction by shocks is not well known 
and we need to resort to a simple parametrization. If, for 
example, the destruction efficiency is very low and nearly 
all the formed dust grains survive, are dispersed and re- 
main in the IGM, the cumulative amount of dust could 
be significant. We choose the dust depletion factor to be 
/dpi=0.4 (Schneider et al. 2003a), and assume a low de- 
struction efficiency, edcst=0.1. Using the global formation 



rate for Population III stars in equation (|13|) and assuming 
/pisn = 1, we calculate a conservative upper limit for the 
mass density of dust in the IGM, 

fi*£M w 10 -7 (25) 

by z = 15. This is much smaller than, for example, the 
density of intergalactic dust, fijjust ~ 10~ 4 ~ 5 at z = 0.5 
(Aguirre 1999; see also Pei, Fall, & Hauser 1995; Croft 
et al. 2000), and thus appears to make negligible contri- 
bution to the dust in the diffuse IGM even if the grains 
survive until the present epoch. We point out, however, 
that even a small quantity of dust in the gas accelerates H2 
production and hence affects the efficiency of gas cooling 
and star formation when the gas metallicity is still low (Hi- 
rashita et al. 2002). Todini & Ferrara (2001) argue that 
H2 production on dust grains dominates over the usual H~ 
channel when the dust-to-gas ratio is greater than 5 per- 
cent of the local Galactic value D g ~ 6x 10~ 3 . Our model 
predicts that the mean dust-to-gas ratio is ~ 3 x 10~ 6 
at z = 15 for / dpl =0.4 and e<iest=0.1, still an order of 
magnitude smaller than the critical value. H2 production 
on dust grains could still be important if the dust-to-gas 
ratio were locally high, either because the distribution is 
inhomogeneous, or because the dust grains remain prefer- 
entially within a confined region around the SN site. 

3.6. Building-up an Early X-ray Background 

It has been suggested that X-rays could act as a source 
of positive feedback by boosting the free electron fraction 
in distant regions. This would in turn promote the forma- 
tion of H2, hence enabling enhanced cooling of primordial 
gas (e.g., Haiman, Rees, & Loeb 1996; Oh 2001; Glover & 
Brand 2003). Apart from early quasars (e.g., Eisenstein 
& Loeb 1995; Bromm & Loeb 2003; Madau et al. 2003), 
SNRs are the only plausible X-ray sources at very high 
redshifts. High redshift SNRs cool via inverse-Compton 
cooling, losing a large fraction of their energy to CMB 
photons (Oh, Cooray, & Kamionkowski 2003). Some frac- 
tion of the thermal energy is also converted into X-rays 
via thermal bremsstrahlung, and yet another unknown 
fraction goes into relativistic electrons which could emit 
X-rays via inverse-Compton scattering of CMB photons. 
Unfortunately, the relative fractions of energy delivered to 
these processes is highly uncertain, and likely to be time- 
dependent during the evolution of the SNRs. In order 
to obtain a rough order-of-magnitude estimate, we cal- 
culate the evolution of an early X-ray background using 
some working assumptions as follows. We consider only 
X-rays with energies greater than 1 keV, because the IGM 
is nearly optically thin to these X-rays, and assume fur- 
ther that a fraction cx of the total SN explosion energy is 
converted to radiation in the relevant X-ray band, between 
1 keV and 5 keV. We set, somewhat arbitrarily, cx = 0.1 
as our fiducial value. Note that this value is significantly 
larger than the one assumed by Glover & Brand (2003) 
as a maximum efficiency case, based on the estimate of 
Terlevich et al. (1992) for X-ray bright remnants in dense 
environments, cx ~ 0.01. Hence, the results presented 
below provide a conservative upper limit. 

For the emissivity of the SNRs, we rather crudely model 
it to be described by a power law L v cx v in the narrow 
(1-5 keV) energy range. By noting that radiative cooling 
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in the SNRs occurs mainly in the first few million years 
for the original SN energy of ~ 10 53 ergs (see section 2.3), 
the time-averaged luminosity in X-rays from a single SNR 
is estimated to be 

L v pa 4.2 x 10 20 1 ergs s _1 Hz _1 , (26) 

where hv§ — IkeV. We compute the evolution of the diffuse 
X-ray flux by solving the cosmological radiative transfer 
equation (e.g., Peebles 1993) 

where J is the specific intensity in units of ergs s _1 cm 
Hz -1 sr _1 , H(z) the Hubble parameter at redshift z, c the 
speed of light, and e the proper volume-averaged emissiv- 
ity We use this unconventional unit for the X-ray flux 
in order to facilitate comparisons with the commonly con- 
sidered soft-UV background radiation. We compute the 
volume emissivity by combining equation (|26|) with our 
model SFR (equation [13], see also Figure 1), assuming 
that all the Population III stars die as pair-instability SNe 
(/pisn = !)■ This simple model again serves as a maximal 
case in terms of X-ray emission efficiency. 

Figure [3] shows the evolution of the X-ray background 
radiation intensity at 1 keV. It is clear that the flux reaches 
only up to 10 -24 ergs s _1 cm~ 2 Hz -1 sr _1 by z — 15. 
Machacek, Bryan, & Abel (2003) performed numerical 
simulations of pre-galactic object formation to examine 
the effect of the positive feedback from a diffuse X-ray 
background in the early universe. They conclude that the 
net effect is quite small for X-ray fluxes below ~ 10 -23 
ergs s _1 cm~ 2 Hz -1 sr -1 . To result in a positive feed- 
back that counteracts the negative feedback due to photo- 
dissociating radiation, an order of magnitude higher flux 
is necessary, which is two orders of magnitude higher than 
our estimate for the X-ray flux at z = 15. Overall, the 
effect of the early X-ray background caused by SNRs to 
the primordial gas cooling appears unimportant at z > 15. 



3.7. The G-dwarf Problem 

It has often been suggested that the paucity of metal- 
poor stars in the solar neighborhood relative to the predic- 
tion of theoretical models, the so-called G-dwarf problem, 
can be resolved by invoking early massive Population III 
stars (Cayrel 1986; see also Larson 1998 for a recent re- 
view). In fact, the kind of pre-enrichment by early gener- 
ation stars proposed here may offer an excellent resolution 
if the metallicity of the gas is sufficiently high initially. 
The simplest models require a significant degree of pre- 
enrichment in the Galactic disk (e.g., Binney & Tremaine 
1987), Zi n it = 0.25^0, whereas our model predicts far too 
small a metallicity accomplished by Population III stars 
alone, ~ 10~ 4 Zq. Thus, considering the self-regulating 
nature of early star formation (see §3.1), we conclude that 
massive Population III stars do little to resolve the classi- 
cal G-dwarf problem. 




redshift 

Fig. 3. — The evolution of the X-ray background radiation inten- 
sity. We use unconventional units, ergs cm 2 s _1 Hz - 1 sr — 1 , in order 
to compare with the typical intensity of the soft-UV background ra- 
diation Juv ~ io -21-23 (in the above units), which is responsible 
for the negative feedback on H2 cooling. 

However, massive Population III stars may account for a 
modern version of this problem - the paucity of extremely 
metal-poor stars in the Galactic halo. As we have shown 
in previous sections, a prompt enrichment to ~ 10 Zq is 
achieved by z — 15. In fact, the gas metallicity around 
the SN sites, where most plausibly the second genera- 
tion stars will be formed, is higher than this mean value, 
and thus, within our model, the second generation stars 
are expected to have metallicities at or above this critical 
value. Observationally, with the important exception of 
the ultra-metal poor star discovered by Christlieb et al. 
(2002), the paucity of stars with [Fe/H] < —4 is apparent 
in our Galaxy. It would be extremely important to derive 
the abundance pattern of such stars through a systematic 
search (e.g., Christlieb 2003). 

3.8. Running Spectral Index Models 

The rapid metal-enrichment by massive Population III 
stars has important implications for some variant CDM 
models. Based on a joint analysis of the first-year WMAP 
data, the 2dF galaxy redshift survey, and the Lyman-a 
forest, Spergel et al. (2003) conclude that cosmological 
models with a running primordial power spectrum are fa- 
vored over power-law models with n « 1 (see also Peiris 
et al. 2003). However, post- WMAP studies indicate that 
the observed high optical depth is not achieved in such 
models with a strong running, because the formation of 
pre-galactic size objects is systematically delayed owing to 
the suppressed power on small scales (e.g., Somerville et 
al. 2003). Chiu, Fan, & Ostriker (2003) and Avelino & 
Liddle (2003) indeed showed that strong running of the 
primordial power results in considerable delay of reioniza- 
tion. Yoshida et al. (2003b) carried out cosmological hy- 
drodynamic simulations of early structure formation and 
concluded that the running spectral index (RSI) model 
is incompatible with the observed high Thomson optical 
depth. However, Yoshida et al. also argue that the result- 
ing Thomson optical depth can be made marginally con- 
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sistent with the WMAP data, if the ionizing photon pro- 
duction rate is extremely high in galaxies at z > 6. This 
practically indicates that, in the RSI model, a top-heavy 
IMF for the stars formed in galaxies must be maintained 
over a long period, 6 < z < 20. 

The problematic nature of this requirement can be ap- 
preciated by noting that the total elapsed time to z = 6 
in a ACDM cosmology is approximately 10 9 years. Very 
massive stars, like those considered here, have lifetimes 
~ 2 x 10 6 years. Thus, if star formation commences 
at high redshifts, z £ 20, nearly 500 generations of these 
high-mass stars could have formed. As we indicate below, 
the metals from these stars must be "hidden" to an im- 
probably high degree of precision for Population III star 
formation not to self-terminate well before z = 6. 

Here, we show that, in scenarios such as the running 
spectral index model, a large quantity of metals is neces- 
sarily produced if it is assumed that Population III star 
formation continues until z — 6. Let us consider a sim- 
plified model in which the stellar IMF is essentially a 5- 
function at 3OOAf , i.e., all the stellar baryons are con- 
verted into very massive stars. This is the most efficient 
way of getting ionizing photons out of a certain mass of 
baryons converted into stars, because the luminosity per 
stellar mass is essentially a constant for metal-free stars 
with mass exceeding 300 Mq. In this case we may simply 
assume a high ionizing photon production rate per stellar 
mass of 1.6 x 1O 48 /M . As in Yoshida et al. (2003b), 
the star-formation rate density for the RSI model is com- 
puted using the same functional fit (equation ^21) with 
p = 1.5, q = 3.05, r = 1.85, a = 0.01, = 0.057, and 
p*(0) = 0.0176M Q yr" 1 Mpc~ 3 . 
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Fig. 4. — We plot the cumulative mass of stellar baryons, which 
is locked up in black holes in our extreme model (top), and the 
mass- weighted mean metallicity (bottom) for the RSI model on the 
assumption that all the stellar baryons are processed through very 
massive stars. 

Figure 0] shows the cumulative mass density of stellar 



baryons. It directly corresponds to the mass density of 
black holes if we assume that all the stars become BHs, 
i.e., their masses are greater than ~ 300M Q . In the plot, 
the dotted line indicates the present-day mass density in 
black holes in our local universe estimated by Merritt & 
Ferrarese (2001). Clearly the mass density increases be- 
yond the present-day mass density by z — 10. On the 
other hand, if all the stars are massive, but within the 
pair- instability SN range of 140M© < M < 26OM , up 
to half of the stellar mass can be processed into heavy 
elements and expelled. Clearly these models are ideal- 
ized, but they allow us to explore the possible parameter 
space. On these assumptions and using the SFR for the 
RSI model, we compute the mass- weighted mean metallic- 
ity. The bottom panel of Figure 0] shows the evolution of 
the mean metallicity Z. Remarkably, the critical metallic- 
ity ~ 10~ 4 Z© is achieved already by z ~ 14 and it reaches 
~ IO^Zq by z — 5. It is now clear that, if we resort to 
the extremely top-heavy IMF in the RSI model to produce 
enough UV photons to reionize the universe early on, the 
total amount of heavy elements quickly rises up to a point 
where the gas chemical properties are expected to change 
(Bromm et al. 2001b). It is conceivable that all the mas- 
sive stars collapse to black holes so that the IGM remains 
chemically pristine without metal-pollution, but then per- 
haps too much baryonic mass goes into black holes, which 
appears to be incompatible with the recent estimate of the 
local mass density in black holes. In summary, a running 
of the primordial power spectrum seems to be incompat- 
ible with various observations in one way or another. A 
similar conclusion is expected to obtain for other variant 
CDM models, such as those involving warm dark matter 
(e.g., Somerville et al. 2003; Yoshida et al. 2003c). 

4. SUMMARY 

We have revisited the cosmological consequences of an 
early generation of Population III stars. We have primar- 
ily considered very massive (~ 100 — 300Af Q ) stars and 
used a detailed model of the global star formation rate 
which is derived from semi-analytic modeling and the re- 
sults of cosmological hydrodynamic simulations that in- 
clude radiative feedback effects. We have shown that the 
first stars alone can produce a sufficient number of pho- 
tons to reionize a significant volume of the IGM by z ~ 15, 
if the stars are very massive, in agreement with cosmo- 
logical radiative transfer calculations (e.g., Sokasian et al. 
2003a). Furthermore, on the assumption that the major- 
ity of the stars explode as pair-instability SNe, the mean 
IGM metallicity increases to the critical value above which 
the gas cooling efficiency is greatly enhanced. This may 
lead to the formation of ordinary stellar populations in- 
cluding low-mass stars (e.g. Mackey et al. 2003). We also 
discussed possible implications of the results. Although 
a broad range of observations indicates that the kind of 
massive stars we considered in the present paper existed 
in the early universe, the simple model may not be able to 
explain the observations in details. Clearly, further studies 
using more sophisticated models are necessary to examine 
the cosmological consequences of early star formation. In 
particular, we will explore the contribution from the sec- 
ond (and subsequent) generation stars forming at z < 13 
in future work. 

The evidence for an early generation of very massive 
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stars may well be strengthened by future observations. 
The ongoing operation of WMAP will yield a more pre- 
cise value for the total optical depth to reionization. In 
the longer term, post- WMAP CMB polarization experi- 
ments such as Planck will probe the reionization history 
(e.g., Kaplinghat et al. 2003). Detection of the second- 
order polarization anisotropies on arcminute scales can 
place strong constraint on the details of reionization (e.g., 
Liu et al. 2001). Near-infrared observations of afterglows 
from high-redshift gamma-ray bursts can also be used to 
probe the reionization history at possibly z > 10 (Barkana 
& Loeb 2003; Inoue et al. 2003; Yoshida & Bloom 2003, 
in preparation). Mapping the morphological evolution of 
reionization may be possible by observations of rcdshifted 
21cm emission by the Square Kilometer Array and LOFAR 
(e.g., Tozzi et al. 2000; Furlanetto, Sokasian, & Hern- 
quist 2004). Data from these future observations will pro- 
vide a more complete picture of cosmic reionization and 
will enable us to distinguish the sources responsible for 
reionization. The precise measurement of the near-IR cos- 
mic background radiation will constrain the total amount 
of light from early generation stars (Santos, Bromm, & 
Kamionkowski 2002; Salvaterra & Ferrara 2003). Ulti- 
mately, direct imaging and spectroscopic observations of 
high redshift star clusters by the James Webb Space Tele- 
scope will probe the evolution of stellar populations up to 
z ~ 10 - 15 (e.g., Stiavelli, Fall, & Panagia 2003; Tumlin- 
son, Shull, & Venkatesan 2003). 



Finally, measurements of the relative abundances of var- 
ious heavy elements in metal-poor stars should provide 
valuable information on the chemical history of the uni- 
verse (e.g., Burris et al. 2000; Norris et al. 2002). Interest- 
ingly, a strong argument against very massive (> 140Mq) 
stars comes from the observed abundance pattern of C- 
rich, extremely Fe-deficient stars (Christlicb et al. 2002; 
Umeda & Nomoto 2003; but see Schneider et al. 2003b). 
It remains to be seen whether or not such stars are truly 
second generation stars and their elemental abundances 
should precisely reflect the metal-yield from the first SNe. 
Observations of a large number of extremely metal-poor 
stars will construct better statistics (Norris, Ryan, & Beers 
2001) and improve constraints on any models for the early 
chemical evolution. Understanding the origin of the first 
heavy elements in the universe and the nature of the 
sources that are responsible for cosmic reionization will 
require the concerted use of data from these broad classes 
of observations. 
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